Temperature Induced Spin Density Wave in Magnetic Doped Topological Insulator 
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We study the magnetic properties of Bi2Se3 doped with isoelectronic magnetic impurities. We 
obtain that at zero temperature the impurities order ferromagnetically, but when raising the tem- 
perature the system undergoes a first order phase transition to a spin density wave phase before 
the system reaches the paramagnetic phase. The origin of this phase is the non-trivial dependence 
of the spin susceptibility on the momentum. We analyze the coupling of the non-uniform magnetic 
phase with the Dirac electronic system that occurs at the surfaces of the topological insulator. 
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I. INTRODUCTION 

Topological insulators (TI) are a newly discovered type 
of systems which are insulating in the bulk and charac- 
terized by the existence of a robust helical gapless Dirac 
two dimensional electron system at their surface [TH3]. 

TI's are typically band insulators for which strong spin 
orbit coupling produces an inversion of the bulk band 
gap. Therefore, in TFs the energy gap is related with the 
spin orbit coupling and that limits its magnitude. The 
most studied and more promising topological insulator is 
Bi2Se3, which is a three dimensional TI with a relatively 
large bulk energy gap '^ 0.3eV and with the Dirac point 
of the surface states located outside the bulk bands [31[S]. 
Angle resolved spectroscopy [HIS] and scanning tunneling 
microscopyfjlexperiments have shown the Dirac nature 
of the surface states of Bi2Se3. 

The spin and wavevector of the surface states of a 
TI are strongly coupled, and the occurrence of a half- 
quantized Hall effect when a energy gap opens at the 
surface has been predicted[8l [9]. Due to the protected 
character of the Dirac states, a gap at the surface should 
be opened with a perturbation that breaks the time re- 
versal symmetry. This can be done by doping the sys- 
tem with magnetic impurities. At the surface of the TI, 
because of the large spin-orbit coupling, the interaction 
between the Dirac-like surface states and the impurities 
induces a large single ion magnetic anisotropy and po- 
larizes the spin of the impurities perpendicularly to the 
surface. This spin-orbit coupling translates in the open- 
ing of an energy gap at the Dirac point of the surface 
states [T0Hl9]. 

From the experimental side, angle resolved photoemis- 
sion spectroscopy (ARPES) studies on the surface of Fe- 
doped Bi2Se3 single crystals have confirmed the opening 
of an energy gap at the Dirac point [20J and the creation 
of odd multiples of Dirac fermions 21 . Also, recently, 
experiments in thin films of Cr-doped Bi2,Sb2-a;Te3, has 
shown a large anomalous Hall conductance in a magnet- 
ically doped topological insulator [22]. 

However, recent experiments |23j found that the spins 
of Fe ions deposited on Bi2Se3 orient inplane. Also 



ARPES experiments [211 ES] found Dirac crossing even 
in the presence of magnetic impurities in contradiction 
with earlier experiments and existing theory. On the 
other hand, recently it has been reported the suppression 
of the Dirac point spectral weight, both in magnetically 
doped and undoped TI, suggesting that the observed gap 
at the Dirac point can not be taken as the sole evidence 
of a magnetic gap [26]. In addition, density functional 
theory based calculations ^2T find that Co adatoms lying 
in the Bi2Se3 surface exhibit an energetically stable mag- 
netic moment perpendicular to the surface, whereas for 
Co atoms located on the interlayer Van der Waals spacing 
the momentum is in the plane parallel to the surface. All 
these results indicate the complexity of the interpretation 
of the ARPES experiments and the possible importance 
of other effects not included in the Dirac hamiltonian, 
as crystalline anisotropy or surface reconstruction, might 
play an important role on the orientation of the magnetic 
impurities. In this work we use an effective hamiltonian 
for describing Bi2Se3, which although it does not include 
microscopic details of the material describes appropri- 
ately the basic properties of the Bi2Se3 related with its 
band structure topology. 

In this work we study the phase diagram of magneti- 
cally doped Bi2Se3. Bi2Se3 is a layered material formed 
by five atom layers arranged along the z-direction. We 
find that at low temperatures the magnetic impurities 
order ferromagnetically along the z-direction. By raising 
the temperature, the TI undergoes two transitions; A 
first order transition from the ferromagnetic to the spin 
density wave phase, and at higher temperatures a sec- 
ond order transition from the spin density wave phase to 
the paramagnetic phase. The spin density wave phase 
has both, the polarization and the wavevector, parallel to 
the z-direction. We have also studied the effect of the 
surface states by calculating the magnetization as func- 
tion of temperature of a slab of Bi2Se3 topological insula- 
tor. Here we find that the surface magnetization survives 
to higher temperatures than the bulk spin density wave 
phase. 

The paper is organized as follow, in Section II we de- 
fine the hamiltonian we use for describing the electri- 
cal properties of Bi2Se3. In Section HI we calculate the 



wavevector dependent paramagnetic spin susceptibility 
of Bi2Se3 and discuss the interaction between magnetic 
impurities through the paramagnetic susceptibihty. In 
Section IV we formulate a Landau theory for describing 
the magnetic order of magnetically doped Bi2Se3, and 
discover the existence of a ferromagnetic to spin density 
wave phase transition at finite temperature. In Section 
V we study the polarization profiles of a magnetically 
doped Bi2Se3 slab and analyze the effect that the Dirac- 
like surface states have on the magnetic phases. We finish 
in Section V with some conclusions and remarks. 



II. HAMILTONIAN 

The low energy and long wavelength electronic proper- 
ties of Bi2Se3 topological insulators are described by the 
four bands k • p HamiltoniaufS], 



H= E{\<i) 



(1) 



+ Al(k)T^ (g) I+Aik^T,^ ® (Tz+A2[k^T.j. + kyTy) ® a,^ 

where an and r^ and Pauli matrices, / the unity 
matrix, Al(k)=Mo - B2{kl + k^) - Bikl, k±^k^ ± 
iky and E{k)=C + Dikl + D2{kl + k^). The 
Hamiltonian is written in the basis |1 >=|plj,t>, 
|2 >=-i|p2-,t>,|3 >=|pl+,i>, |4 >= i|p2j,i>, which 
are the hybridized states of the Se-p orbital and the 
Bi-p orbital with even (+) and odd (— ) parities and 
spin up (t) and down (J,). The Hamiltonian param- 
eters for Bi2Se3 are[28] , Mo=0.28el^, Ax=Q.22eVnm, 
A2=0AleVnm, Bi=0.10eVnm'^, B2=0M6eVmn'^, C=- 
0.0068ey, Di=0m3eVnm^ and D2=0.196eVnm^. In 
this basis the spin operators get the formpQ]. Sz — 
I ® Uz, Sx — Tz® Ox and Sy — Tz® Oy. 



III. BULK SPIN SUSCEPTIBILITY 

The paramagnetic susceptibility obtained from the 
Hamiltonian EqfTlhas the form 



x/xM(q) 
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E E 



n occ. 
n empty 



|<n',k + q|5'^|n,k>| 



(2) 



Here \n, k > and e„,k are the eigenfunctions and eigenval- 
ues of Hamiltonian EqJTJand Jl is the sample volume. In 
the case of an insulator, this spin susceptibility is caused 
by the coupling of the valence and conduction band in- 
duced by the spin operator pOI- The susceptibility is a 
smooth function of the wavevector and because the sys- 
tem is an insulator there are no anomalies associated with 
Fermi surfaces. The symmetry of the original Hamilto- 
nian dictates that the non-diagonal elements of the sus- 
ceptibility tensor are zero and Xxx = Xyy ¥" Xzz- 

In FigjTlwe plot the Xxx and Xzz as a function of qz and 
qx- The direct coupling, A2k±, between atomic orbitals 
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FIG. 1: (Color online) Spin susceptibility as function of 
the wave vector along the z and the x-directions. The arrow 
indicates the position of the maximum. 



with opposite parities and opposite z-component of the 
spin, makes that for k± ^ 0, occupied and empty states 
are coupled through Sz- Whereas those states are only 
connected through Sx when kz 7^ 0. This makes x^z(q) > 

Xxx{.(\). 

The more important contribution to Xzziqz) comes 
from regions in the reciprocal space where the matrix 
elements <n' jk + qz\Sz\n,'k> |, with n occupied and n' 
empty, reaches the maximum value. This happens when 
A^(k)=0 or MCk + qz)=0- For a given kz these con- 



ditions define two circular crowns of radius 



B~2 



and 



Mo-Bi{k,+q,)^ 



B2 



and thickness v42/(2i?2)- There- 



fore, the area of the reciprocal space that contributes 
appreciably to Xzz{qz) increases with qz- For larger val- 
ues of qz one of the circular crowns collapses to zero and 
the contributions to the integral decrease. This behavior 
explains qualitatively the maximum that Xzz presents at 
a wavevector G ^ ^/Mq/Bi. 

The existence of a maximum in Xzzi^lz.) at finite qz 
is robust against small changes in the parameters of the 
four bands Hamiltonian. In FigJ2]we plot Xzz(qz) for dif- 
ferent values of the TI gap. The position of the maximum 
decreases continuously towards q — when Mq decreases 
and only disappears for small values of Mq. In the normal 
insulator phase, Mq < 0, the maximum always occurs at 
g = 0. 



A. Coupling between diluted magnetic impurities 

Consider now a TI doped with magnetic impurities of 
spin iS*. We assume that the number of electrons in the 
system does not change in the presence of the magnetic 
impurities. That can be achieved by doping with iso- 
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crystal approximation [321 [331 IMl I3S|, and we consider 
that the system is invariant in the {x,y) -plane, and the 
polarization only depends on the z-direction. In the next 
section we obtain the magnetic polarization as function 
of temperature and z-coordinate by using a Landau free 
energy functional. 



IV. LANDAU FREE ENERGY FUNCTIONAL 

We assume that the system is invariant in the (x, y)- 
plane, and allows the polarization to oscillate with period 
2tt/G along the z-direction. In consequence we define the 
normalized magnetic polarization, —1 < m[z,T) < 1, as, 



FIG. 2: (Color online)Spin magnetic susceptibility, Xzz, for 
different values of AIq, as function of the momentum in the 
z-direction. AH other parameters of the band structure cor- 
responds to those of Bi2Se3. As the value of the mass Mq 
decreases the position of the maximum of the susceptibility 
moves towards small values of q^ ■ At small values of Mq the 
maximum occurs at gz = 0. In the normal insulator phase. 
Mo < 0, the maximum of x^z occurs at q^ — for all values 
of the mass parameter. The dots indicate the position of the 
maximum. 



electronic magnetic dopants or by adding compensating 
non-magnetic dopants [3T]. In this work we consider the 
dilute limit, i.e. concentration of impurities smaller than 
5%, for which the direct interaction between the spin of 
the magnetic impurities can be neglected. 

However, the electrons spins have a strong exchange 
coupling, ^J(|r|), with the magnetic impurities spins 
which, in turn, are equally affected by the exchange 
field of the electrons. In this form the magnetic impuri- 
ties in the system interact mediated by electronic states. 
We treat this interaction in second order perturbation 
theory [321 [33, that has proved to be a reliable approx- 
imation in diluted magnetic semiconductors [331 134|. In 
this approach the effective exchange parameter between 
two magnetic impurities separated by a vector R and 
spins pointing in the z^-direction is. 



J.(R) 



S^ 



j2^/J7^X..(q)e^''^ 



(3) 



where Jeff = / J(|r|)dr is the effective exchange cou- 
pling between the magnetic impurity and the electron 
spin. 

Because Xzz > Xxx in all range of wavevectors, the 
system has an easy axis of magnetization along the z- 
direction and therefore isoelectronic magnetic impurities 
in Bi2Se3 will tend to polarize in the z-direction. The 
maximum that the spin susceptibility presents at finite 
wavevector in the z-direction, will determine the exis- 
tence of non-uniform polarization in magnetically doped 
TI. We treat the magnetically ordered state in the virtual 



m{z, T) = moiT) + mciT) cos(Gz) 



(4) 



where niQ and thq are the order parameters of the uni- 
form ferromagnetic (FM) phase and the spin density wave 
phase (SDW) respectively. 

The internal energy per unit volume corresponding to 
this magnetization is, 

E = -'^mlx.ziO) - '^mlxzziG) ■ (5) 

where J = ^J^ffC, being c the density of magnetic im- 
purities. In our case, the value of Xzz{G) is less than 
10% larger than Xzzi^) and the zero temperature ground 
state is a uniform FM phase, m{z,T = 0) = 1. However 
the maximum of the spin susceptibility at G will modify 
the spin density at larger temperatures. 

Knowing that for small values of the polarization, the 
entropy of a classical spin at a given T is, see Appendix 
A, 

-TS^-kBT\n{2) + ^kBTm^ + ^kBTm^ (6) 

we get that in the mean field approximation and for small 
values of magnetic polarization, the Landau free energy 
per unit volume takes the form. 



J^=-2'^oXz 



:(0) - -^m^GXzz 



(G) 



1 1 



\n2 -^m^{z,T)-^m\z,T). 



(7) 



where /3 — l/ksT and L is the sample dimension in the 
z-direction. Using the expression Eq[4| and in the limit 
L — >■ 00, we get 



J" = -mlksiT - To) 



-m^kBiT -Tg) 
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(8) 



where Tq and Tq are the critical temperatures of the pure 
FM and SDW phases respectively. 



kBTo'- 



J 



[xzziO) and ksTo^^XzziG) . (9) 



The phase diagram of a system described by a free energy 
as that of Eq|8] depends on the relative magnitudes of the 
fourth-order potentials [H^ . In our case the product of the 
pre-factors of ttiq and m^ is smaller than the square of 
the rn^rriQ pre-factor and there is no phase coexistence 
in the phase diagram. By increasing the temperature, 
there is a first order transition from the FM phase to the 
SDW phase at 



rri* 



^/3To - V2Tg 



(10) 



This is the main result of this work: by heating, a mag- 
netically doped TI undergoes two phase transitions, a 
FM to SDW first order transition at T* and a SDW to 
paramagnetic second order transition at Tq. Although 
at T=0, the FM phase has lower energy than the SDW 
phase, the FM to SDW transition at finite T occurs be- 
cause the entropy of the SDW increases faster with T 
than the entropy of the FM phase. 

In the next section we analyze how the surface states 
existing in topological insulators couple to the bulk mag- 
netic polarization. 



V. SPIN POLARIZATION OF MAGNETICALLY 
DOPED TI SLABS 

At the surface of a TI there exists a two dimensional 
Dirac electron gas. Because the chirality of the electron 
gas, an exchange field perpendicular to the surface opens 
a gap in the spectra. Then, in order to minimize the en- 
ergy, a magnetic impurity will polarize perpendicularly 
to the surface P^UHT^ [T5] . In the diluted limit, surface 
states mediate an RKKY interaction among the impuri- 
ties which is always ferromagnetic, whenever the chemi- 
cal potential resides near the Dirac point pUl W]\ [TSl 155] . 
Therefore magnetic impurities at the surface of a TI will 
order ferromagnetically perpendicular to the surface. 

We are going to study numerically the spin polarization 
as function of temperature and position of a magnetically 
doped TI slab. The objective here is first to confirm 
the results obtained with the Landau functional where 
we consider a unique Fourier component of Xzz{(lz) ^nd 
second to analyze the coupling between the surface and 
the bulk magnetization. 

We analyze a TI slab of thickness L and perpendicular 
to the z-direction. We expect the electron affinity of 
Bi2Se3 to be much larger than its band gap. Therefore, at 
the surface of the TI we will neglect the penetration of the 
electron wavefunction into the vacuum. The eigenvalues, 
Era.kj and wavefunctions, ^n,k(-2^)i are obtained by solving 
Eqll] with kz — ~idz and forcing the wavefunction to 
vanish at z = and z = L. This is satisfied expanding 
^n,\i{z) in harmonics. 
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FIG. 3: (Color online)In the left panel we plot the band 
structure of a topological insulator film 30nm thick. For this 
thickness there is no coupling between localized states on op- 
posite surfaces and surface states are degenerated. Dashed 
areas denote the bulk band structure. In the right panel we 
plot the absolute value of the four components of the wave- 
function of a surface state with momentum close to zero. 



here A is the sample area and we choose N,nax large 
enough so that the results does not depend on it. 

For L > lOnm, the surfaces of the slab are decou- 
pled and the band structure is independent of L. In the 
bulk energy gap region, appear some surface states which 
are the benchmark of the TI. In Figj3] we plot the band 
dispersion and the shape of the wavefunction of these 
states. The results we obtain agree completely with pre- 
vious results [31 [29] . 

In the slab geometry the momentum in the z-direction 
is not a good quantum number and the paramagnetic 
susceptibility depends on two position indices z and z' . 
Therefore, in the virtual crystal approximation and in 
second order perturbation theory, the internal energy of 
the magnetically doped TI slab is. 



E- 
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2L 



L pL 



dzdz'x{z, z')m{z)m{z'), with (12) 



X{z.z) = - J2 






K.Mi^^)Sz^n'Mi^^) X *;, ,^(z')5.*„,k(/) , (13) 

where npi^) is the Fermi distribution function. xi^T^') 
indicates the coupling between uniform polarized (x, y)- 
planes, located at positions z and z'. The interaction 
between magnetic impurities is mediated by electrons in 
the system, and because the bulk system is an insulator, 
the interaction is very short ranged in the z-direction, see 
Figg 

We compute the temperature dependence of the mag- 
netization profile in the mean field approximation. At 
a given position z, the magnetization m{z,T) feels a (in 
energy units) magnetic field, 



B{z) = J dz'xiz,z')m{z') , (14) 

Jo 

and the magnetization of an isolated impurity in the pres- 
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FIG. 4: (Color onlinejXzz {z,z') evaluated at the maximum 
of the surface wavefunction, z=1.2nm, and at the center of a 
20nm thick slab, as function of z' . The first case corresponds 
to a region near the surface, where the two dimensional Dirac 
electron system contributes to the response functions. In the 
latter case the response function is not affected by the sur- 
face and it is the bulk response function. In both cases the 
functions are very peaked at z=z' . The negative values of 
the coupling in the bulk response function is a consequence of 
the maximum that the response function present at qz=G in 
the reciprocal space. Near the surfaces, and because of their 
metallic character, the magnetic coupling is stronger. This is 
reflected in the asymmetry of the dashed line, the interaction 
between planes is larger as closer the planes are to the surface. 
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FIG. 5: (Color online)Magnetization as a function of temper- 
ature for TI slabs of thickness (a) L = 5nm, (b) L = lOnm and 
(c) L = 30nm. To is the bulk FM critical temperature of the 
topological insulator. The small "step" in the middle of the 
first order transition that occurs at L = 30nm is consequence 
of an interference effect between the surface magnetization 
and the bulk SDW phase. 



FIG. 6: (Color online) Magnetic polarizations as a function 
of the position across the topological insulator slab z, for dif- 
ferent layer thickness and temperatures. 



ence of the molecular field is. 



m{z, T) = coth 



B{z) 



B 



(15) 



Solving self-consistently Eq[T4] and EqjT5J we obtain the 
magnetization profiles as a function of T. 

Because the metallic surface states intermediate a 
RKKY coupling [ini El] at the surface, the response func- 
tion, x{z, z') is larger near the surface than in the bulk, 
see Fig|4j Therefore, as function of T, the absolute value 
of the magnetization decreases faster in the bulk region 
than in the surface [5^. However, it is important to note 
that the surface and the bulk are part of a unique sys- 
tem and therefore there is only a unique critical tempera- 
ture, corresponding to the transition of the paramagnetic 
phase. 

In FigjSlwe show the magnetization as function of tem- 
perature for TI slabs of thickness L = 5nm, L = lOnm 
and L = 30nm. We plot the average value of m{z,T), 
and the value of the magnetization on top of the surface 
states. In Figure l6J we plot the magnetization profiles 
for different temperatures and L = 5nm, L = lOnm and 
L = 30nm. 

For L ~ lOnm and L =^ 30nm the surfaces are practi- 
cally decoupled and the central part of the slab behaves 
as bulk. There is a strong jump in the magnetization at 



T*, that indicates the first order FM to SDW transition. 
In the SDW phase the oscillating magnetization does not 
contribute to the total magnetization and the magnetiza- 
tion for T > r* is due to surface states. In Fig|6Fb)-(c), 
it is apparent at the center of the slab, the abrupt tran- 
sition from an uniform magnetization phase to a SDW 
phase. For smaller thickness of the slab, FigJ6Fa), the 
surface states are coupled and there is no well defined 
bulk region, that reflects in the absence of FM to SDW 
transition. 

The magnetization at the surface is practically not af- 
fected by the FM to SDW transition, and decays with T 
continuously to zero. The ferromagnetism at the surface 
is more robust than in the central part. For tempera- 
tures where mp and uiq are near zero, the surface of the 
system can be more than 30% polarized. These results 
indicate the possibility that the magnetization at the sur- 
faces of TFs could be finite at temperatures larger than 
the bulk critical temperatures Tq and T* |39j . Because of 
the metallic character of the TI surface states, there is a 
range of temperatures, for which the Dirac-like electron 
system at the surface of the TI is gapped, although the 
bulk part of the system is practically unpolarized. 

A similar SDW phase has been also obtained numer- 
ically by Rosenberg and Franz in a slab geometry of 
Bi2Se3[32]- However these authors interpret the oscil- 
lation of the polarization as spatial fluctuations of the 
bulk magnetization coupled with the surface magnetiza- 
tion. From our calculation we attribute the oscillations 
in the magnetization reported in ref . [3 9) as a signature 
of the bulk SDW phase. 



VI. FINAL REMARKS AND CONCLUSIONS 

In this work we study the phase diagram of magnet- 
ically doped Bi2Se3. At low temperatures the magnetic 
impurities order ferromagnetically along the z-direction. 
By raising the temperature, the TI undergoes two tran- 
sitions; A first order transition from the ferromagnetic 
to the spin density wave phase, and at higher tempera- 
tures a second order transition from the spin density wave 
phase to the paramagnetic phase. This results could ex- 
plain recent experimental results |40j that suggest the ex- 
istence, as function of the temperature, of two different 
magnetic phases in Fe doped Bi2Se3. 

We have also studied the effect of the surface states by 
calculating the magnetization as function of temperature 
of a slab of Bi2Se3 topological insulator. Here we find 
that the surface magnetization survives to higher tem- 
peratures than the bulk spin density wave phase. The 
existence of a range of temperatures for which the bulk 
magnetization practically vanishes whereas a finite mag- 
netization exits at the surface, could explain some ex- 
perimental results that observe a gap at the surface of 
Bi2Se3 but not bulk magnetism [20) [21] . 

It is important to analyze the behavior of the phase 
diagram as function of the gap parameter Mq. In FigjT] 
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FIG. 7: (Color online) Phase diagram of a thick magnetic 
doped TI, as function of the mass parameter Mq. 



we show the phase diagram of a magnetically doped thick 
TI slab as a function of Mq. For Mq < the system 
is a normal insulator and there are no surface states. 
Also the spin orbit coupling is small and the SDW phase 
does not exist. For Mq > the system is a TI and 
the gap increases with Mq. TI with larger gaps have 
more metallic surface states and the FM order at the 
surface is therefore more robust. Also the effective spin 
orbit coupling is stronger and both T* and Tq increase 
with Mq . The results of FigJT] show that the range of 
temperatures where the SDW phase exists increases with 

Mq. 

Finally we make an estimation of the critical temper- 
ature. From the band structure parameters of Bi2Se3, 
choosing the density of the magnetic impurities to be 
5xl0^°cm~'^, the total angular momentum of a single 



magnetic ion to be S" = 3/2 and the effective exchange 
coupling Jeff =250meVnm^^2\ we obtain r^""= k18K. 
This values can change by factors of two by changing 
the magnetic ions or the density of impurities. It is well 
known that the mean-field approximations tend to over- 
estimate the transition temperature due to the neglect 
of the fiuctuations. In diluted magnetic semiconductors 
thermal fluctuations reduce the value of the Curie tem- 
perature in near 30% [33], and we expect a similar reduc- 
tion in topological insulators. 
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Appendix A: Mean Field Expression for the Entropy. 

In this section we obtain an expression for the entropy 
of a system of classical spins of magnetization m, that are 
coupled with the topological insulator trough a general 
term _E[to]. 

The Free Energy of a system of classical spins of mag- 
nitude unity in an external magnetic field h is 



1 /smM^X 
(3 \ Ph J 



(Al) 



from where the magnetization can be calculated as 



m =< m >= 



9J- 



1 



1 



dh tanh(/3/i) /3/i ' 
The entropy of the spin system is then 



(A2) 



(A3) 



The total energy of the system is 

1 



J" 



total 



E[m]-TS = E[m]- 



smr 



(A4) 
where E[m\ is the change in the electronic energy of the 
system because of the polarization of the magnetic im- 
purities. 

To obtain h, we minimize the total free energy with 
respect to h, dJ^°^^^ /dh = 0. In the hmit of small h 



sinh(^ (^_(^ 

^ jih ' ^ ' & 180 



(A5) 



In this limit, the minimization condition gives ph 
— 3to — ym? -\- ..., and the entropy gets the form 



TS=-kBT\n{2) + ^kBTm^ + —k^Tm'^. (A6) 
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